Directed Sub- Wavelength Imaging Using a Layered Metal-Dielectric System 
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We examine some of the optical properties of a metamaterial consisting of thin layers of alter- 
nating metal and dielectric. We can model this material as a homogeneous effective medium with 
anisotropic dielectric permittivity. When the components of this permittivity have different signs, 
the behavior of the system becomes very interesting: the normally evanescent parts of a P-polarized 
incident field are now transmitted, and there is a preferred direction of propagation. 

We show that a slab of this material can form an image with sub- wavelength details, at a position 
which depends on the frequency of light used. The quality of the image is affected by absorption 
and by the finite width of the layers; we go beyond the effective medium approximation to predict 
how thin the layers need to be in order to obtain subwavelength resolution. 

I. INTRODUCTION 

An anisotropic material in which one of the components of the dielectric permittivity tensor has a different sign 
to the others has interesting properties. It supports the propagation of modes which would normally be evanescent, 
and these modes travel in a preferred direction. The propagation of evanescent modes gives us hope that an image 
produced by light travelling through a slab of such a material might retain a sharp profile; also, because the preferred 
direction depends on the ratio of the components of the permittivity tensor, it can be controlled by varying the 
frequency of light used. 

We first look at a way of producing a metamaterial with the desired properties: by making a system of thin, 
alternating metal and dielectric layers. A system of this type was proposed by Ramakrishna et al.i as a form of 
"superlens" ; it improves on the original suggestion for a superlens^ which consists of just a single layer of metal, and 
has recently been realisedi&i 

We then look at the dispersion relation for our anisotropic material, to see why modes which would be evanescent in 
both the metal and the dielectric separately are able to propagate in the combined system, and why there is a preferred 
direction of propagation. The subwavelength details of the source are transmitted through the system because they 
couple to the surface plasmons^ that exist on the boundaries between metal and dielectric; this mechanism is the 
basis for the current interest in metallic structures for super-resolution imaging at optical frequenciesi^iiiiiSil 

Next, we investigate the transmission properties of a slab of this material, and apply our formulae to the case of a 
line source. We show that we can expect to obtain a sharp image as long as the amount of absorption is not too high. 

Finally, we go beyond the effective medium approximation to show the effect of the finite layer widths on the optical 
properties. We demonstrate that the "resolution" of the slab is limited by the width of the layers; thinner sheets 
mean that the description of the system using the effective medium becomes increasingly accurate, and the image 
quality improves. 

II. LAYERED SYSTEMS 



We concentrate on periodic layered systems of the form shown in figure ^ We assume that each layer can be 
described by homogeneous and isotropic permittivity and permeability parameters. When the layers are sufficiently 
thin, we can treat the whole system as a single anisotropic medium with the dielectric permittivit y 9 1 10 
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FIG. 1: System geometry. The layers are infinite in extent in the iy-plane. 



A helpful way to see this is through the characteristic matrix formalism^ this method, which is related to that used 
by Rytov in the original derivation^ is described in the appendix. 

The homogenized magnetic permeability is given by expressions analogous to and When rj is small, the 
effective parameters arc dominated by the first medium, while for large 77, they resemble those of the second medium. 

Only the ratio of the thicknesses of the two layers appears in the homogenized version, not the absolute value; 
however, the characterization of the material using the effective medium parameters is more accurate when both d\ 
and g?2 are small. 

For a layered metal-dielectric system, we can tune the response either by altering the frequency or by changing 
the ratio of layer thicknesses. 11 This is demonstrated by figures [21 and which show the real and imaginary parts of 
the effective permittivity for two different thickness ratios, for a system composed of alternating layers of silver and 
silica. The material data from which these plots are constructed have been taken from the books by Palik— and 
Nikogosyan— . In both graphs, there are two regions in which Si(e x ) and SR(e z ) take opposite signs. In the first region, 
which includes energies up to approximately 3.2 eV, 3t(e x ) is negative; in the second, which consists of a small range 
of energies around 3.6 eV, 3?(e x ) is positive. 

By choosing a suitable value of 77, we can make the real parts of e z and e x take opposite signs over a range of 
frequencies. We investigate the consequences of this in the next section. 



The unusual behavior of the layered materials can be understood by considering the dispersion relation between 
the frequency to and the wave vector k. We assume that we are dealing with non-magnetic materials, so that the 
magnetic permeability /i = 1. If the dielectric permittivity is anisotropic, the interesting waves are those with 
transverse magnetic (TM) polarization. The dispersion relation for these waves is 



We have taken k y to be zero, since the x- and y-directions are equivalent. When e x and e z are both positive, the 
relationship between k x and k z is similar to that in free space: for small k x , k z is real, but when k x becomes large, 
k z becomes imaginary. The propagation of the wave in the z-direction is governed by k z ; when k z is imaginary, the 
wave is evanescent: it decays exponentially with z. 

However, when e x and e z have opposite signs, k z is real for a much wider range of values of k x . Even the high 
spatial frequency components with large k x , which would normally be evanescent, now correspond to real values of 
k z , and hence to propagating waves. 

If we want to plot the dispersion relation, we have to remember that the permittivity itself is frequency-dependent. 
To get an idea of what the dispersion relation looks like, we can use an idealized model: we imagine a metamaterial 
whose layers are composed of equal thicknesses of a dielectric, with positive, frequency-independent permittivity, and 



III. PERMITTIVITY WITH DIRECTION-DEPENDENT SIGN 
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FIG. 2: (Color online) The dielectric permittivity of the metamaterial constructed from layers of silver and silica. This and 
the following graph show the real and imaginary parts of the in-plane and perpendicular components of the permittivity for 
different layer thickness ratios; in this case, 77 = 2/3, which means that the layers of silica are one and a half times as thick as 
the layers of metal. 



a metal, with the simple plasma-like permittivity 

e m (w) = e m (oo) j. (5) 

For now, we assume that the materials are non-absorbing. The resulting dispersion relation is plotted in figure We 
can identify two distinct bands from the figure. In the lower, e x is negative, while e z is positive; the signs are reversed 
in the upper band. In both cases, the contours of constant to are hyperbolae. In the lower band, these hyperbolae are 
centered on the k x -axis, while in the upper, they are centered on the k z -axis. 

In fact, there is also a third band at high frequencies, but this is the least interesting regime and is not shown in 
figure both components of the permittivity are positive here. 

The dispersion relation also provides the key to the preferred propagation direction. This is determined by the 
group velocity. A constant- frequency section of the dispersion relation (cut across the first band) is plotted in figure 
[5j The hyperbolic form of the curve means that for large \k x \, it tends to the following straight line: 




The group velocity is perpendicular to the constant-w contours like the one plotted in figure|3] The figure demonstrates 
that apart from a region around k z = 0, the group velocity vectors all point in almost the same two directions: this 
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FIG. 3: (Color online) The effective permittivity when 77 = 1.5. The silver layers are now thicker, and the amount of absorption 
has increased: the imaginary parts of the permittivity are now larger. However, the real parts are also correspondingly larger 
in magnitude. 



is the basis for the preferred direction of propagation. Remembering that the x- and y-directions are equivalent, we 
can see that the preferred directions form a cone around the z-axis. The half-angle of the cone is 



9 — arctan. /— — . (7) 
V £z 

In the region around k z = 0, the arrows point outside the cone. In this band, there are no propagating modes in a 
small region around k x = 0, and no propagating modes with a group velocity vector lying inside the cone. 

If we take a cross-section from the second band, instead of the first, we also see a hyperbolic contour; the plot 
resembles figure [5] but rotated by 90°. The group velocities for the modes around k x = now point inside the cone, 
rather than outside. 

To conclude this section, we look at the physical process that allows our layered metamaterial to mimic an anisotropic 
material and to support the propagation of normally evanescent waves. The key fact is that surface plasmons are 
supported at an interface where the permittivity changes sign. When the metal permittivity is negative, this sign 
change occurs at every interface; the wave is transmitted via coupled surface plasmons, as indicated in figure 
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FIG. 4: (Color online) The dispersion relation for an idealized metal-insulator system. The permittivity of the metal is given 
by |J5J with e m (oo) = 2.0, while the dielectric has permittivity td = 2.5; the layers are of equal width (rj — 1). k p is the wave 
vector corresponding to the plasma frequency (fc p = ui p /c). The first two bands are shown; they have been separated to make 
visualisation easier, but there is no band gap. The plots are symmetric about the planes k x = and k z — 0. 

IV. TRANSMISSION THROUGH AN ANISOTROPIC SYSTEM 

We have seen that we can produce a metamaterial with interesting properties by stacking alternating layers of metal 
and dielectric. Next, we look at a slab of this material, and examine the transmission coefficient. 

We assume that the slab is embedded in a uniform medium of constant permittivity (which may be unity, repre- 
senting vacuum). In such a medium, the dispersion relation J2J becomes 

k 2 x + kf = k 2 e. (8) 

We write k' z to distinguish the z-component of the wave vector in the surrounding medium from that in the slab. The 
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FIG. 5: (Color online) The relationship between k x and k z for lu = 0.2ui p (in the middle of the lower band in figure QJ. The 
straight lines show the asymptotes given by equation ©. The group velocity is indicated by the arrows, which are perpendicular 
to the curve; the length of the arrows is proportional to the magnitude of the group velocity. 




FIG. 6: (Color online) Schematic diagram of the transmission of normally evanescent waves, showing the role of surface 
plasmons. The line represents the electric field strength. 
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FIG. 7: (Color online) The transmission coefficient, as denned in equation with the surrounding medium taken to be air. 
The slab width used is d = 2/k p , which would be of the order of lOOnm for a plasma frequency cu p ~ 4eV. The results for three 
different frequencies are plotted, corresponding to the three bands referred to in the discussion of figure |1] The real part of e x 
is negative when u) = 0.3ui p , while that of e z is positive; the signs are reversed when lu — 0.65cj p . At higher frequencies, both 
are positive. The permittivity of the metal is taken to be e m = 1.7 + 0.6i — Up/aJ 2 , while that of the dielectric is ed = 2.5. 



transmission coefficient for TM waves is 

Kk x , w) = (9) 

2 cos k z d - i [^f + |t J sin k * d 

where the dispersion relations Q and iJSJl are used to define k z and k' z in terms of k x and to. 

In figure \7\ we plot the transmission coefficient for three different regimes, corresponding to the three different 
frequency ranges. At high frequencies (here represented by u> — oj p ), both components of the metamaterial permittivity 
are positive. In this regime, the transmission coefficient is close to unity for small wave vectors. It drops abruptly 
to zero at k x — ko, and rises equally sharply afterward, again approaching unity; finally, it decays exponentially for 
larger wave vectors. Very similar behavior is observed in the intermediate frequency range (lj — 0.65oj p ). This time, 
the maximum following the zero at k x = ko is higher, and the rate of exponential decay for large k x is less rapid. 

The most interesting frequency range is the lowest one (w = 0.3w p ). There is the usual zero in the transmission at 
k x = ko, followed by a very sharp peak. However, there is also significant transmission even for large wave vectors; 
the transmission coefficient has a series of peaks, decreasing in magnitude, and approximately periodic in k x . The 
resonances correspond to localized states for the slab; they are in turn antisymmetric and symmetric. There is a 
difference between the first two resonances (just above k x = ko) and those for higher wave vectors. For the first two, 
the wave is non-propagating inside the slab (because k z is almost purely imaginary): the resonances therefore consist 
of coupled surface plasmons located on each surface of the slab. For the higher wave vectors, the wave is able to 
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FIG. 8: (Color online) Imaging a solenoidal line source. 

propagate^ within the slab (because k z is almost purely real), and the transmission peaks correspond to Fabry-Perot 
resonances - standing waves inside the slabii 

In fact, a similar set of peaks would be visible in the intermediate-frequency regime, were it not for absorption. 
The material parameters used to generate figure \J\ include a realistic amount of absorption, and a glance at figures [21 
and |3 shows that absorption is high in the region where 9?(e z ) becomes negative. The localized states are supported 
in both the low- and intermediate-frequency ranges, but are suppressed in the latter by high absorption. 

V. IMAGING A LINE SOURCE 



We have seen that the layered system allows enhanced transmission of high-spatial-frequency components at certain 
frequencies. This gives us hope that we may achieve sub- wavelength imaging using the slab As a test, we consider 
the image of the line source pictured in figure |S| 

In the absence of the metamaterial, the field generated by this source is 



E(r)= f [ e lk * x+zk y y+lk '* z - wt E a (k y )(it-^z) dk x dk y , 



(10) 



where the current profile in the y-direction is as yet unspecified. As before, k' z represents the ^-component of the 
wave vector in the surrounding medium. 

When we place the metamaterial next to the source, as shown in the figure, some radiation will be reflected from the 
slab and will generate additional currents. If we neglect these, we can estimate the ^-component of the transmitted 
field as 



E™(r) = ^ r e^+iWKiz-D-^H (Ji% + k*,w) ■ E °j ky ]\ x dk x dk y 

J -oo J -oo ^ V ' k x + k 



(11) 



Note that this is the TM component of the field. In general, there will also be a TE component which must be 
calculated. However, if we consider a line source which is uniform in strength and infinitely long, so that E (k y ) <x 
S(k y ), the entire field is transverse magnetic. In this case, the calculation reduces to the solution of the following 
integral: 



El M (r) 



E 



k o J-oo 2 cos k z d - i {^^- + -j^f-^j sin k z 



(12) 
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FIG. 9: (Color online) The transmitted field intensity, comparing analytical and numerical results. The parameters used were 
£d = 2.5, e m = 1.7 + 0.05i — ujp/oj 2 , with U — 0.68u p . In each case, we plot the field just beyond the slab (so that for example 
when k d = 2.0, k z = 2.01). 



The integral can be solved approximately when the frequency is in the intermediate range: that is, when 3ft(e x ) > 
and 3?(e z ) < 0. The resonant states then all have large k x ; they are the standing wave states discussed in the previous 
section, rather than the coupled surface plasmon states (which have k x close to ko). We are therefore justified in 
making the near field approximation, which leads to the following analytic form for the ^-component of the transmitted 
field: 

TM ±4niE e -ifci(| a |-«(»-d))-«^ 

x [V) ~ k Q d(l/e z -e x ) ' i + e -iAfc(M-i(*-d)) ' 



where 

fc- = 

d V e x V 1 + e x(z 



k x = -> arctan ( h/IJ^k I , 1 4 ) 



and 



In this approximation, i?J M and Ej M are identical to within a phase factor. In figure we plot the intensity of 
the transmitted field, comparing the approximate analytical solution with the results of numerical integration. To 
generate the plot, we take an unrealistically low value for the absorption in the metal; the point of the graph is to 
compare numerical and analytical results, but also to demonstrate the features which we hope to be able to observe. 
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FIG. 10: (Color online) A schematic showing that reflections lead to periodically-repeated images of the two principal peaks. 

First, we note that the position of the peaks is proportional to the slab width. This is a manifestation of the 
preferred direction of propagation; within the metamaterial, the light travels at a fixed angle to the z-axis, in the 
xz- plane (since we have translational invariance in the y-dircction). The secondary peaks which are visible when 
kod =1.0 are caused by reflection from the boundaries; this is why they overlap precisely with the primary peaks for 
the slab with k d — 3.0. The reflections are illustrated in figure ITUl 

In the first frequency regime, the approximate analytical solution is more difficult to obtain: there are the additional 
surface plasmon resonances close to k x = k$, for which one cannot make the near held approximation. However, it is 
still possible to obtain numerical results. As one would expect from figure [71 these are much more promising: using 
realistic parameters, we are able to produce a sharp image, as shown by the line marked "Effective medium" in hgure 
1111 The width of the principal peak in the effective medium approximation is around A/10. Figure ITT1 also illustrates 
the results of a more detailed analysis, which goes beyond the simplified effective medium approach; we will discuss 
these next. 



VI. BEYOND THE EFFECTIVE MEDIUM APPROXIMATION 

Treating the layered system as an effective medium is a helpful simplification, in terms of both understanding and 
performing simulations. However, it has limitations. In this section, we model the system in more detail, considering 
the hnite width of the layers; naturally, as the layers are made thinner, we see that the effective medium approximation 
becomes more appropriate. 

First, we look at the dispersion relation for the layered metamaterial. We can obtain k z as a function of k x (at a 
given frequency) from the characteristic matrix, as described in appendix B These isofrequency contours are plotted in 
hgure IT^l The effective medium has full translational symmetry, but this is broken when considering the structure of 
finite-width layers; the system becomes periodic, and the figure shows part of the first Brillouin zone, which extends 
from k z — — Trfcp to k z = irk p . The effective-medium contours are deformed by the new periodicity, and bend towards 
the zone boundaries. This introduces a new cutoff: for a given frequency, there is a value of k x above which no 
propagating solution exists. This affects the resolution of the lens-like system. 

The next logical step is to investigate the change in the behavior of the slab of metamaterial described in section 
IIVI From now on, we focus on the first band. Figure 1131 shows that the new cutoff in k x is clearly manifested 
in the transmission function: above the cutoff, the transmission decays very rapidly. Below the cutoff, we see the 
familiar Fabry-Perot and coupled surface plasmon resonances, although they have moved slightly; this is because the 
relationship between k x and k z has been altered, as shown in figure 1721 

Finally, we re-examine the image of a line source using the modified transmission functions shown in figure 1131 
Figure 1771 shows the transmitted electric field intensity, plotted as a function of x, for various different layer widths. 
Increasing the width of the layers which make up the metamaterial slab (while keeping the total slab width constant) 
causes the principal peak to broaden, as expected. As the layers get thinner, the transmitted image more closely 
resembles the effective medium result. 



VII. CONCLUSION 

We have investigated a class of anisotropic materials in which the one of the components of the dielectric permittivity 
has a different sign from the others. These materials are able to support the propagation of modes that would normally 
be evanescent: they are able to collect and transfer the near field. In addition, inside the anisotropic medium, light 
travels in a preferred direction. 
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FIG. 11: (Color online) The transmitted electric field intensity for a line source, imaged by a metamaterial slab of thickness 
l/fco- The material parameters used correspond to layers of Ag and ZnS-SiCh, embedded in crystalline Ge2Sb2Te5 (a phase- 
change material used in optical storage devices), for light of wavelength 650nm. This corresponds to a total slab width of 
around 105nm. 



We have studied the transmission properties of a slab made up of such a material. The image of a line source 
consists of two lines, with an offset determined by the ratio of the components of the permittivity; the width of the 
imaged lines depends on the amount of absorption, but in principle can be much less than the wavelength of light 
used. 

One realization of such a material is a stack of alternating layers of metal and dielectric. The thinner the layers, the 
better this metamaterial approaches the form of the ideal anisotropic medium. We have shown how the ideal band 
structure is deformed by the non-zero layer width. Using realistic material parameters, we have also demonstrated 
that a stack of alternating Ag and ZnS-SiC>2 layers can form an image of a line source which is much narrower than 
the wavelength of light when working at 650nm. 

The combination of subwavelength resolution with the fact that the position of the image depends on the frequency 
of light being used suggests that this layered system may have useful applications. For example, in conjunction with a 
super-resolution near- field optical structure (supcr-RENS)^*i^ it may allow the possibility of multiplexed recording. 

APPENDIX: HOMOGENIZATION IN LAYERED SYSTEMS 

The effective medium parameters for our one-dimensional system of alternating layers can be calculated by using 
the characteristic matrix method. The effective dielectric permittivity is obtained from a consideration of TM fields; 
TE fields give the effective magnetic permeability. 
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FIG. 12: (Color online) Isofrequency contours demonstrating the effect of finite layer width. The solid lines are the contours 
for a system where the cell size (di + d,2) is l/fc p ; the result is a Brillouin zone of extent 2-k on the fc z -axis. The dashed lines 
are the equivalent contours in the effective medium approximation. The three sets of contours correspond to the three bands 
of the dispersion relation discussed in section UTTl the first two are hyperbolic in the effective medium approximation, while the 
third is elliptical. The material parameters are the same as those used to generate figure [I] 



The geometry of the system is shown in figure ^ with the z-axis perpendicular to the layers. We take the plane 
of incidence to be the irz-plane; the symmetry of the system means that this is equivalent to the yz-plane, and the 
results which follow are general. 

The characteristic matrix Mj (k x , d) relates the Fourier component of the field in the plane z = zq to that in the 
plane z — zq + d (all within medium j). For TM waves in a homogeneous medium, it takes the form£ 



Mj(k x ,d) 



cos ki 3 d 



sin&a d coskz'd 



i k c j 

IF 



sin ki^d 



(A.l) 



where kl is given by the dispersion relation 



(A.2) 



The matrix for a single cell of our layered system, consisting of one sheet of each material, is just the product of 
the matrices for the separate layers: 



M cea (k x ,di,rj) = M 1 (k x ,d 1 )M 2 (k x ,r)d 1 



(A.3) 
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FIG. 13: (Color online) The transmission as a function of k x for various layer widths. The total slab width is maintained at 
1/fco in each case, while the number of individual layers is adjusted. 



A stack of n cells has the characteristic matrix M„ = (M co n) n . We can calculate this by diagonalizing M ce n; we then 
obtain 

We have introduced A, which is one of the eigenvalues of M cc u; the other eigenvalue is A -1 , which follows because 
det M cc n = 1. We have also introduced p and q, which are the ratios of the components of the eigenvectors of M co u. 

Expanding in powers of the layer thickness allows us to relate M„ to the characteristic matrix for an effective 
medium: 



M « = U^-- ^ ^ " )+0[{l + v)dx] (A.5) 



coskfn(l + ?7)c?i ^Lsinfc^ ff n(l +rj)di 
sin fc| ff n(l + r))di cos k° s n(l + if}d\ 
= M cS + 0[(l + r 1 )d 1 ] (A.6) 
where the effective medium parameters are 

ei + V e 2 



1 + V 

1 _ 1 rj 

V z ~ £1(1 + 7;) + 62(1 + 7?)' 



(A.7) 
(A.8) 
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These parameters appear in the dispersion relation for the effective medium, which differs slightly from 1A.2|) because 
the permittivity is now anisotropic: 

^ + ^J_ = k l (A.9) 

We also note here that the cell matrix M CG n has another use. We can determine the true dispersion relation for the 
layered system - without using the effective medium approximation - by finding the eigenvalues and eigenvectors of 
this matrix. When the eigenvalue has unit modulus, we have found a Bloch mode; we then make the association 

A = e ife,(i+^)di_ ( A10 ) 
The eigenvalue A depends on the frequency and on k x ; equation (|A.10|I is therefore the dispersion relation. 
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